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Abstract 



We propose the difference discrete variational principle in discrete mechanics 
and symplectic algorithm with variable step-length of time in finite duration 
based upon a noncommutative differential calculus established in this paper. 
This approach keeps both symplicticity and energy conservation discretely. 
We show that there exists the discrete version of the Euler-Lagrange coho- 
mology in these discrete systems. We also discuss the solution existence in 
finite time- length and its site density in continuous limit, and apply our ap- 
proach to the pendulum with periodic perturbation. The numerical results 
are satisfactory. 
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I. INTRODUCTION 



In the course of numerical calculations in continuous mechanics, the equations of motion 
should be discretized in time variable to get certain schemes. Among various discrete time 
schemes in Hamiltonian mechanics, the symplectic algorithm is quite successful [15,4,16], 
since it preserves the symplecticity. In Lagrange mechanics, the discretizations may be 
established via discrete variation method to get the symplectic-preserving integrators. But 
this variation is usually taken with respect to the discrete dependent variables only [10], [18]. 

Thus, for the long time, there had been no discrete Legendre transformation. Therefore, 
the discrete Lagrange and Hamilton mechanics could not be transformed to each other. In 
addition, the equal-time-step symplectic integrators [18,14] cannot preserve the energy, while 
the variable time-step variational integrator proposed in [10,11] preserves discretely energy 
but without proving the symplicticity. Until recently, it is proposed [9,1,2] that the energy 
conservation and symplicticity for the variational integrators with variable time-steps may 
be compatible. However, it is still open how to deal with discrete Hamilton systems as well 
as symplectic algorithm by variational principle. In the meantime, the difference discrete 
variational principle (DDVP) has been introduced in [6,7]. And in this approach the discrete 
Legendre transformation can be introduced and both discrete Lagrange systems including 
variational integrators and Hamilton systems including the schemes in symplectic algorithm 
can be transferred to each other. One of the key points in DDVP is that the difference is 
treated as an independent variational variable. Very recently, DDVP with variable time step 
(VDDVP) has been proposed in [5]. Thus the symplicticity and energy could be discretely 
preserved at equal footing. 

In this paper, we propose the non-commutative differential calculus (NCDC) on fibre 
bundle over discrete base space in order to deal with the variation of some Lagrange system 
with high order derivatives. Based upon this NCDC we present the VDDVP with variable 
time steps as well as the discrete version of the Euler-Lagrange cohomology and apply them 
to discrete mechanics including the generalized symplectic algorithm. For a given discrete 
action S within finite time-duration, we obtain both discrete equations of motion and energy 
conservation equation from the total variation 5S := i^ total d t S = 0, where d t is the total 
exterior differential operator, by means of VDDVP with fixed end points of the duration. 
Then we show that the discrete version of the first Euler-Lagrange cohomology does exist in 
the case of variable time-step. We also consider some simple example to show the advantages 
of our approach. 

The paper is organized as follows. In section 2, we propose an NCDC on the fibre bundle 
whose base space is 1-dimensional lattice with variable step-lengths. In section 3, we present 
the VDDVP and deal with both discrete Lagrangian and Hamiltonian mechanics based upon 
such an NCDC. In section 4, we discuss the solution existence in finite time step-length and 
site density in continuous limit. In section 5, by means of the VDDVP, we obtain energy- 
conserving midpoint scheme for Hamiltonian system that may depend on time manifestly 
and show its numerical properties. Finally, we end with some remarks and discussions. 

II. THE NCDC ON LATTICE WITH VARIABLE SPACING 

For the lattice with equal spacing, an NCDC on its function space has been proposed in 
[8]. In this section, we generalize it to the one on the lattice with variable spacing in order 
to describe the VDDVP on discrete mechanics including energy-preserving and symplectic 
algorithm. 
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A. The NCDC on Lattice with Variable Spacing 

For the simplicity, we consider one dimensional case. It corresponds to the time variable 
t G R and is dicretized as t k G To, 

* g R^t g T D = {(t k ,t k+1 = t k + At k , he I)}. (2.1) 

Note that At k are not fixed here. The total number of nodes in set / is finite for the case 
of finite interval, while it is is infinite if the intervalis infinite. We will discuss both of 
them. The algebra of all the functions defined on Tp is denoted as A. There is an algebra 
homomorphism defined by the shift operator as follows 

R E :A^A, R E f(t k ) = f(t k+1 ). (2.2) 

Then one introduces the vector fields V = span{dE} on A, where 

(d E f)(t k ) ■= (Re - id)f{t k ) = f(t k+1 ) - f(t k ). (2.3) 

The space of 1-form f2 1 (T D ) = span{x E } is dual to V: x E (9 E ) = 1- The whole differential 
algebra fi*(T D ) can be defined as tt*(T D ) Vt n (T D ) with A = fl°(T D ). Introduce the 

n=0,l 

exterior differentiation in fi* d : fi°(T D ) -> Q}(T D ). It acts on a 0-form / G fi°(T D ) = A 
gives 

df := (d E f)x E e O 1 . (2-4) 

An NCDC with respect to d on A is constructed [8]. As a simplified version of [3], its 
nilpotency and Leibniz law of its acting on forms lead to the noncommutative property in 
X E f(a) = (R E f(a)) X E . 

It is convenient to consider a set of auxiliary bases of Q 1 : 

dt k := (t k+1 - t k ) x E - (2.5) 

Thus the exterior differentiation (2.4) becomes df(-) = [A k f(-)] dt k , where / G A and 
A k f{-):={t k+1 -t k )-\R E -id) f{-). 

For each interval [t k ,t k+ i] on lattice To, the discrete integral over it can be defined as 
follows: At moment t k , if there has a 1-form A k f(t k ) dt k , whose integral over the path 
tk —> t k +i is 



/ A k f(t k ) dt k ■= / Vt k A k f(t k )f(t k+1 ) - f(t k ), (2.6) 

JD,t k JD,t k 

where / denotes discrete integral, and Vt k is the measure at moment t k . Under this 
Jd 

definition, any integral over the path U — > tj, i < j, can be divided into (j — i) parts: 

r am = [ tk+i [^ k f(t k )}dt k 

j d,u i<k<j-i JD ' tk (2 7) 

= E (/(**+!) - /(<*)) = - 

i<k<j-l 

Here discrete integral (2.6) can be regarded as the analogues of the ordinary integral in 
continuous case, and the equation (2.7) is the discrete version of Stokes' formula. 
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B. The NCDC on Fibre Bundle over Discrete Space 

In order to establish the differential calculus on fibre bundle over discrete space, we 
consider the base space, time t, is discretized as 

t g R^t g T D = {{t k ,t k+1 = t k + At k , kel)}, (2.8) 

but its n-dimensional fibre space M k , at moment t k {k G I), is still continuous and smooth 
enough. Let N be the set of all nodes on To with index set Ind(N) = I, and M = \J ke j M k 
be the fibre space on Tp. At the moment t k , M k denotes the set of nodes neighboring to t k , 
and I k is the index set of nodes of M k including t k . The coordinates of M k are denoted by 
q l (t k ) = q l ( k \i = 1, • • • , n. Let M k = Uiei k M be the union of fibre spaces M h I G I k on Af k , 
and F(M k ) be the function space. 

The exterior differential calculus on F(Ai k ) consists of two parts: the vertical differen- 
tiation d v along the fibre spaces and the horizontal differentiation dh along the base space 
[8]. Note that the operator d v acts on fibre variables only, while the operator dh also acts 
on function space F(T D ). Obviously, on the function space F(T D ), the operator dh should 
just play the role of operator d, which is already established in [8]. Thus, at moment t k , on 
total space M. k x To, the total differential operator on function space F(A4 k x To) locally 
is d t := d v + dh- When the functional Lagrangian under consideration depends on time t 
manifestly, we should consider that dh acts not only on F(A4 k ), as the discrete differen- 
tiation with respect to t k via the variables on fibre spaces A4 k , but also on F(T D ) as the 
differentiation directly for variables t k . 

For the local coordinates q l on the fibre, there are 

d t q^ : = d v qW+d h q^ k \ , , 

dnqW := (R E q m - q m ) X E = (q i(k+1) ~ X* \ 1 J 

where Re is translation operator on the fibre bundle space. For the / G F(Ai k x To), the 
operator Rf acting only on fibre spaces is also needed. They are defined as follows, 

Re /({?}; {t k }) ■= f({R E q};{REt k }), 
R f f({q};{t k }) ■-= f({R E q};{t k })- 

The symbol {/} denotes the set of variables of same kind with the function /. For example, 
{q} is the set of variables on fibre space, and {t k } is the set of variables on base space, etc. 
Then we can define the differential calculus on fibre bundle over discrete base space, 
Definition 2.1: On the total space M k x T D , the NCDC on function f G F(M k x T D ) 
is given by : 

dtf({q}; {h}) := E ( V>/) ^ + (R E f- Rf f) X E - (2.11) 

i,ke{q} 

Now, let us construct the higher order differential algebra. As is well known, d v is an 
ordinary differential operator satisfying the ordinary Leibniz rule and acting on F(M) only. 
Its differential complex is also as usual and denoted as fi*(M)0fi n (M). The differential 

n 

complex for d t is the topological product of Q*(M) and fi(To). The differential complex 
of Q(Td) is noncommutative and with only one generator \ E ■ Any functions and Q,(Td) is 
not commutative. Acting d v on the noncommutative relation \ E Q 1 ^ = q l( - k+1 ^ x E ■> there is 
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(2.10) 



~X E A d v q 1 ^ = d v q l ( k+1 ^ A \ E ■ It can be generalized to any differential algebra uo on the 
fibre bundle: 

X E A lu = (-l) de9 ^ (R E cu) A X E ■ (2.12) 
The operation of total operator d t on d v q is defined as follows 

d v d v q i<yk ^ := 0, 

d h d v qW ■= -(R E d v qM - d v qW) A X E ■ 

It is obvious to get that d v \ E = and dh = since there is only one 1-form on base space 
and it is antisymmetric. We can get, 

d v d h qM = (d v q^ - d v qW) A X E (9 ^ 

= -4 d v q^. {2 - U} 

Thus, for the local coordinates q % of fibre, a well-defined differential calculus has been 
established and the total differential operator d t = d v + dh is indeed nilpotent: 



d 2 t = (d v + d h ) 2 = 0. 



■i _ i ,] ,1 xi _ ^ (2-15) 



C. The NCDC on Extended Bundle over Discrete Space 

In [18], [14], Veselov uses Q x Q for the discrete version of the tangent bundle of a configu- 
ration space Q. The Lagrangian L(q, q) is discretized to L(q^ k \ g( fc+1 )). In his approach, how 
to introduce the discrete Legendre transformation and discrete canonical momentum had not 
been considered. In addition, the time step-lengths are fixed with equal step-length. Very 
recently two of us with their collaborators in [6] , [7] , [5] proposed another discrete variational 
principle and introduced discrete Legendre transformation and discrete canonical momen- 
tum. In terms of NCDC on the time lattice (the base space) with variable step-lengths, 
exterior differentiation rules of variables t k 's and q 1 ^ 's are taken on different ways. We 
will follow this line to discuss the discrete variational principle with high order derivatives 
or differences. 

The continuous Lagrangian L(q\ q\ ■ ■ ■ , (J7) V; t) can be discretized in such a form : = 

L{vf\v?\---,vf> k) -t k ) where vf ] := q^,vf k) := A k q^, „•<*> := A fc tfi,2 < m < I. 
In order to deal with the variation problem of the Lagrangian with I — th order derivatives, 
we first consider that the all variables (vq^ , , • • • , vf^) are independent from each other 
and denote the fibre bundle over one dimensional discrete space as E(M,Te>) , then use 
the constrains to fix the relation between them Vq^ := q l ( h \ and vffl := A fe i^L\, 1 < m < I. 

At each point of t k G T D the fibre space M k is of dimension n x (I + 1) and with local 
coordinates (v$ k \ Vi k \ ■ ■ ■ , vf^). In order to define the exterior derivatives we should discuss 
the functions defined on the nodes neighboring to t k . Similar to the discussion in the last 
subsection we can define M. k = \J P ei k M p be the union of fibre spaces M p , p G I k on Af k , and 
F(M. k ) be the function space. 

Now we start from the case of / = 1, where Lagrangian function depends on first order 
derivative only. In terms of definition 2.1, the total exterior differential of A k q 1 ^ is 
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dt(A k qM) = A fc (^«) + (A fe+1 g l ( fc+1 )-A fc g^+ 1 ))x B 

= A k (d v qW) + (A k+1 q^ - A k q^) X E (2.16) 
= d v (A k qW) + [(R E - id) A k q^} X E . 

This relation suggests that the exterior differential for v should be as follows 

d t v\ ik) := d v v\ ik) + d h v\ (k \ , 
d h v^ := \{R E - id) vf\ X E = ivf +1) - v^)x E . { V 

Here the action of translation operators R E and Rf of (2.10) are extended to the function 
with the new coordinate v \^ is : 

Re fiM; {t k }) := f({R E v 1 }; {R E t k }), R f /({^}; {t k }) : = f({R E Vl }; {t k }). (2.18) 

where R E v\^ := Vi k+1 \ 

We can define the total exterior differential d t on the function which depends on first 
order derivative only as: 

Definition 2.2: the NCDC on function f G F(M k ) can be defined : 

dtf({q}, M; {h}) = E E W) d t vW + (R E f- R f f) X E . (2.19) 

m=0 i,k£{v m } 

Similarly we can define NCDC on the case of function depending on higher order deriva- 
tives as follows. 

Definition 2.3: If f e F(M k ) the NCDC on function f is: 

n 

dtf({q}, M, • • • , K}; {tk}) = E E dA k) + (R E f- Rf f) X E , (2.20) 

m=0i,ke{v m } m 

where 

d t vW : = d v v^+d h vf\ 

d h vW := [{R E - id) v^] X E = (4" +1) - v^) X E (9 91 x 

R E f({vm},{h}) ■= f({R E v m },{R E t k }), {Z - ZL) 
R f f({v m },{tk}) ■= f({R E v m },{t k }). 

Then this operator d t on general functions case is defined, its nilpotency on the differential 
complex can be easily constructed as in the last subsection. 



III. DIFFERENCE DISCRETE VARIATIONAL PRINCIPLE IN DISCRETE 
MECHANICS AND SYMPLECTIC ALGORITHM 

In this section, we will study the variation problems in discrete mechanics including 
the symplectic algorithm completely based on the established NCDC on the extended fibre 
bundle where the functions may depend on high order derivatives. The method we will 
use is the closest analogue of geometric-variational approach to continuous mechanics. One 
of the key point is that the difference as the discrete derivative in the NCDC is taken as 
independent variable in the variation. Therefore, the variation of this kind is called the 
difference discrete variational principle [7,5]. 
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A. Difference Variation in Discrete Mechanics 

As were shown in [5], total variation 5 t S in continuous mechanics can be expressed as 
a contraction between total exterior differentiation of action S and total variational vector 
field C,totai- Variational principle says that it should be zero with fixed end points. On the 
other hand, at the path which satisfies Euler-Lagrange equation, the only contribution to 
StS comes from the endpoint: 5tS = Gfa) — G(t\). This is the most general idea of action 
principle in continuous mechanics. It can easily be generalized to the discrete mechanics. 

When we discuss the variation problem in discrete mechanics using the extended fibre 
bundle with local bundle coordinates (vq k \ v^, ■ ■ • , vf®), it is necessary to take vffl = 
Afcfm-i> m— 1, • • • , I, and q 1 ^ = v l Q k \ But the exterior differentiation of f(vffl) is different 
from that of f(AkV l ^ k l l ) by the NCDC, expect for simple linear function. One way to solve 
this problem is to introduce Lagrange multipliers used at constrained dynamics, then we can 
take the discrete action as: 

N-l t i 

Sd:=J2 W? i( *\«f\...,«f ) ;*0+EE^ fc HA fc t;a-4 fc) )}dt fc , (3-1) 

fe=l JD ' t k m =l i 

here A^'s are the Lagrange multipliers with dimension (n x I) and discrete integral J D is 
defined as in (2.6). At the moment t k , as the exterior differentiation rules for each X 1 ^ is 
taken the same form for vffi, as in definition 2.3 except that dimensions of each extended 
fibre space are doubled. 

The variational vector field at moment tk is defined as follows: 

dl := 6t k d tk +± E (rf ) -Aur, + E E T^y, (3-2) 

m=0i,fc'e{-} OVm m=l i,k'e{-} 

where d tk g : = (tk+i — tk)^ 1 (RE9 — Rfd)- The variations 5t k , S t v 1 ^ and<5 t A^ fc '^ can be defined 
from the contraction between total exterior differentiation of dt k , dtv 1 ^" 1 and dfX 1 ^ and the 
basis of total variational vector field £ tota i as follows: 

St k '■= (tk+l - t k ) < X E , Ztoial > = < dt k , £total >> 

5 tV f) := < d t vf'\ £L >= 6 v v^ + (A k ,vW)5t k ,, m = 0, . . . , n, (3.3) 
5 t XP := < d t \f'\ 41 >= 6 v X$p + (Av>&W)6tv, m = 1, . . . ,n. 

Now, the variation of discrete action comes from two parts, one is from contraction 
between d t {L^ + J2 ^ % m\^-kVm-i~ v nt^)} an d4L> another is from variation of the measure. 
At moment tk, there is S(J^^ Vt k ) = S(Vt k ). In comparison with equation 

5( f tk+1 Vt k ) = 5(t k+1 - t k ) = 8t k+1 - 5t k , (3.4) 

JD,t k 

it means 

5(Vt k ) :=< d (Vt k ), tilt > (Vt k ) A k 5t k . (3.5) 
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B. Difference Variation in Lagrangian Mechanics 

In discrete Lagrangian mechanics, a discrete version of variational principle can be ex- 
pressed as follows. At finite interval [ti,t N ] (finite notes) of discrete time coordinate, the 
only contribution to total variation of discrete action should come from end point term. 

Case 1: the /-order Lagrangian = L(q t( - k \ v . . . , v]^; t k ): 
The total exterior differential of S$ is, 

d t S^ = [ tk+1 {(L^ +EE A^A^e^ - vf^)) d(Vt k ) 

JDt k m=l - (3.6) 

+Vt k d t (L^ + ELi Af )(A fc tfi - vf)))}. 
After a straightforward computation we can get, 

dt {A^(A^S)} = (A^S) + X^Mdtvfl,) - A^+ 1 )(A fe+1 4- + i 1) )A,^ fc , 

dt {^ l i k) vf ] )} = vf) d t \f) + ^d t vf\ 

(3.7) 

Then the total variation of becomes 



' D 
I 



StS^ = f k+1 Vt k {A k ( £ E Ajf - £ (fc) ***) + (A fc £« + 4^ (fc) )«* 

jD *fe m=l i 

+ E( E (A fe 4- } i - + (^ (fc) L( fe ) - A fe Af - 1} ) <W W (3.8) 

j m=l 

+ E (^(^ (fe) - A fc A^i 1} " & fa™ + (^)L (fe) - Af )) ))}, 

rn=l 

where := £Ui(£i AjWtfjW) - L«. 

Variational principle says the bulk part of total variation must be zero. After some 
algebraic computation we finally obtain, 

l—m 

t#*> = Afctffi, Af) = (E(-A fe ^y9 i(fc) )L<*>, m=l,...,/, 

(3.9) 

d qi(k) L^ + (E(-A feJ R^ 1 ) /l 9 i(fc) ) L<*> = 0, A fe i?( fc ) + 4 fc £ (fe) = 0, 



If the time step-lengths are fixed, the last equation of (3.9) has no solution in general 
even if the Lagrangian does not depend on time manifestly. 

Case 2: Lagrangian = L(aq i{k ^ + pq^ k+1 \ v[ {k) ; t k ): 

Let dqiL/W be the derivative w.r.t. variable aq 1 ^ + Pq l( - k+1 ^ in L^ k \ and d q iL^ k ^ be the 
derivative w.r.t. variable aq l( - k ~^ + (3q 1 ^ in L^ k ~ l \ etc. Denote r k := t k+ i — t k . By means 
of relation (3d q iL^5 t q^ k+ ^ = A^Tfc-i/^L^- 1 ^^) + T k 1 T k - 1 (3d qi L( k -V5 t q i ( k \ the total 
variation of S$ becomes 



S t S^ = 5 t f tk+1 Vt k {L^ + Y,^ k \^ m -^ k) )} 

JDt k ■ 

= [ tk+1 Vt k {A k Q% + {AkE (k) + d tk L^)5t k ( 3 - 10 ) 



+ E t ((A k q^ - vl (k) )5 t \W + [L qm ] 5 t qW + (9 l(fe) LW - X^)8 t vf } )}, 
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where 



[L qi (k)} 
E (k) 



= ad q ,m + r k \ k ^ (3d qi L^ - A k (X^), 

= E t X {k) v\ {k) -K k \ (3.11) 

= Ei(X i{k - 1] + r k ^(3d qi L^) <W« - EW6t k . 



According to the discrete variational principle and the definition of discrete integral (2.6), 
the discrete Euler-Lagrange equations are, 

vf ] = A k q* k \ A^jiW, 

[L qi(k) }=0, A k E^ + d tk L^ = 0. 

Let pf^ := c^L^" 1 ) + T k _if3d q iL < " k ~ 1 \ so that the variation equations in (3.12) can be 
rewritten as follows: 

A k p ( t ) = (a + P)d qi L^. (3.13) 
Substitute <9^ into pf +V) = <9„»L (fe) + r k (3d q iL^ k \ we obtain 

° Jk+1) + ^— Q pf ) = d v ,L^. (3.14) 



It is just the discrete Legendre transformation, so that discrete Hamiltonian at interval 
[t k ,t k+1 ] is naturally defined as Jf<*> := E ^'H^V i^ " LiK 

When L( fe ) = |(4 (fc) ) 2 - ^(ag*W + /3g^ +1 )) and a + /3 = 1, it is easy to verify that 
equations (3.13) and (3.14) are just the discrete canonical equations in Hamiltonian system 
H = |p 2 — V(q). There are some interested cases: (i): a — 0, (3 — 1, (ii): a — 1,(3 — 0, and 
(iii):ct = P = \. Both case (i) and (ii) are the known simple symplectic algorithms of order 
1 accuracy, and case (iii) is the midpoint algorithms of order 2. 

Now, let's discuss constraint on a and f3. For simplicity let = |(i>i^) 2 — V(aq 1 ^ + 
f3q l ( k+1 )) ) and r k — > 0, k E Z, then we expand discrete Euler-Lagrange equations in q % W 
compare with Euler-Lagrange equation in continuous mechanics, we have to take a + (3 = 1. 



C. Difference Variation in Hamiltonian Mechanics 

The total variation problem in the discrete Hamiltonian mechanics with variable time 
step-lengths can be dealt with from the viewpoint of NCDC too, where we should express 
the discrete Lagrangian in terms of the discrete Hamiltonian via Legendre transformation. 

The discrete version of variational principle is also expressed as follows: at interval [t±, t^] 
of discrete time coordinate the total variation of discrete action comes from end point term 
only. In Hamiltonian mechanics, if discrete Lagrangian is written as L = pv — H(p,q;t), 
then the discrete action can be taken as J D T>t {pv — H(p, q; t) + A(Ag — v )}. 

Case 3: = Ei(«pf +1) + (3pf ] )v\ {k) - H (k \ where a + (3 = 1 and := H(apf +1) + 
(3pf\aqW + /3g i(fc+1) ;t fc ): 

Introduce Lagrange multipliers A^ fc ) as in case 2. Let denote derivative w.r.t. vari- 
able ap\ k+1 ^ + (3pf^ in H^ k \ and denote derivative w.r.t. variable aq 1 ^ + (3q l{ - k+1 ^ in 
jj(k) reS p ec tively. Its total variation of action at interval [t 1: t k ] reads 
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StSP = f k+1 Vt k {A fe eg } + (A fc 4*> - (fe) )<^ 

JDt k 

+ Ei((ap? +1) + Pp? ] - \^)5 t vf ] + (A k q^ - wjW)^*) (3 15) 

+ Wvf } - HW) + ar fc "Vi(«i H) - ^)) <rf } 



where 



eg 5 := Ei(A i(fe_1) - (3T k - 1 H ( £~ 1 ' ) )8 t q i{ ® 
4* := EiA^^-LW. 

Thus we obtain 

vf ] = A k q* k \ = apf +1) + (3.17) 

and discrete equations of motion, 



p {v m _ H W) + ar -i Tfe _ l( ^-i) _ ^-D) = o, 

A^*- 1 ) + PrrW-iH^ + <*H™ = 0, 



(3.18) 



and the equation for the variable time step-lengths, 

A k E^ - 4 fc ifW = 0. (3.19) 

Which are nothing but the exact symplectic energy-conserving algorithms. 

In terms of relation (3. 16), (3. 18) there is E$ = H^ k \ Moreover, the discrete equations 
of motion (3.18) become 

ar k {A kP f > + If £>) + /^(A^?" 1 ' + = 0. (3 ' 20) 

It should be satisfied at all nodes on lattice T D , so that we have to take 

A k q^ _ H W = o, A k pf ] + Hf = 0. (3.21) 

In this way, we have discrete Legendre transformation — = that transfer 

discrete Hamiltonian to discrete Lagrangian on fibre bundle. And the end point terms 
become standard form <d^ = J2i Pi^^tQ 1 ^ — H^5t k . This case can be regarded as the 
counterpart (in Hamiltonian formulism) of case 2. 

/,From 5 v Sd part in (3.15), it is easy to see that if we only take d v on So, it also can be 
divided into two parts: 

d v S D = £ [ tk+1 Vt k {A k 9 D W + £ D {k) }, (3.22) 
k jD ' tk 

here S D {k \ 6 D (k) are the discrete Euler-Lagrange 1-form and symplectic potential 1-form 
respectively. Due to the nilpotency of d v , it is straightforward to get 

d v £ D ™ + A k uj D ^ = 0. (3.23) 
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Here differential structure : d v 6o^ is just the symplectic structure Y,id v Pi A d v q^ k \ 

Therefore, we may get the discrete version for the theorem of Euler-Lagrange cohomology 
and symplectic structure-preserving condition [6], [7]: 

Theorem: For a discrete Hamiltonian system with H(p,q;t), we have: 

1. The discrete version of the first Euler-Lagrange cohomology may be nontrivial: 

Hdq M := {Closed Euler-Lagrange forms}/ {Exact Euler-Lagrange forms}. 

2. The necessary and sufficient condition for conservation of the discrete symplectic 
2-form, i.e. 

AW° = 0, (3.24) 
is the corresponding discrete Euler-Lagrange 1-form being closed. 



IV. SOLUTION EXISTENCE AND SITE DENSITY 

In the discrete evolution problem, if there has solution whose discrete time series satisfy 
Tfc+i = Tk + 8r k = Tk + 0(r k ), we say there has a stable solution. In this section, we discuss 
the case of H{p, q; t) = |p 2 + V(q, t) with midpoint algorithms, and find this algorithm will 
limit the form of V(q,t) in order that discrete system have stable solution. 

Now, let us analyze the discrete evolution problem of t k — > t k+1 and — > t k+2 . In 
these two steps, we have four discrete canonical equations, 

p( fc+1 ) - p^ _ v , ^ k+r > + qW t fe+ i + t fc g( fc+1 ) - gW _ p( fc+1 > + p^ 

tk+i — tk 2 2 tk+i — tk 2 

p (k+2) _ p (fc+l) q (k+2) + g (k + l) tk+2+tk+1 q (k+2) _ q (k+l) p (k+2) + p (k+l) 



th-LO — t 



= -vc- — ^ — , + , 



k+2 — tfc+l * * tfc+2 — 



(4.1) 



(4.2) 



and two discrete energy equations, one of which is 

1 | t k+1 +t k 

1 p (fc + 2) 2 +p(fc+ l) ( fc+ 2 2) + g(fc+ l) 2 tfc+1+t 

2 l 2 J 1 2 2 J ' 

Above five equations are enough to discuss the relation of r k and 5r k as follows. Firstly, 
let's denote 

t k +i — t k = r k , t k +2 — tk+1 = T k + 5r k , 

p(k+l) _ p(k) _ ^ p(k+2) _ p(k+l) _ r _|_ ^ (4-3) 

_ = ^ g (^ + 2) _ g(fc+l) = s + tfs ' 

In this way, the quantities on the fc-th and (k + 2)-th nodes can be expressed by those 
quantities on the (k + l)-th node and 5r k , 5r and 8s. At the same time, we also use notation 
p _ p(k+i)^ g _ g(k+i) anc j j. _ f or simplicity, so those five equations become 

= -V(q- -,t- -), — = (P- o)> 

T k 11 T k 2 

r + 5r , s + 5s r k + 5r k s + 5s , r + 5r 

-V (q + — , t + ), - ■ - =(p+ — ), 



T k + 5r k 2 2 r k + 5T k 2 

J ' 2 

(4.4) 



!/ r N2 . T/ / S + T ^ l ( . r + ^N2 , T/ / . S + ^ S + T ^ 
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Secondly, If r k is a small time step- length and there are r, s ~ T k and 5r, 5s, ^r^ ~ r|, we can 
expand above equations in Tailor Series and omit its higher order terms than r|. The new 
five equations can be regarded as constraint equations of six variables Tk, r, s, 5rk, Sr and 
5s, so we can finally obtain an equation of Tk and 5T k : 

Vlp r k + \{VV + Wfr + V' 2 ) 5r k + \{Vl'p 2 + \V V) ^ 2 = 0. (4.5) 

If the perturbation term gives us a small parameter ~ Vt/V, which is less than or com- 
parable with a finite small Tk, it should have a solution in the domain V"p 2 + 3V t 'p + V' 2 ^ 0, 
as showed in the next section. In the other words, hypothesis 5rk ~ r| should be modified 
to 5rk ~ s(tk)Tk, where e(tk) is a small quantity but not tend to zero in continuous limit. 
But in the case of finite small Tk, we can keep the form of (4.5) since e{tk)Tk is comparable 
with the 2-order term t| or its higher order term in Tailor series. 

If Hamiltonian does not depend on time manifestly, the hypothesis 5T k ~ r| can be kept 
even in the continuous limit. Now, (4.5) becomes 

( 3VV + 30^ + V-V^0, or ^ = - (3v ^'f 3v , 2) r, (4.6) 
It means 5r k /r k — > too, so we have 

^ = n = n (i + — ) « n e ^ = eE ^° ( 4 - 7 ) 

r fc i=0 T ^+* i=0 Tfc +* i=0 

According to equation (4.6) and definition of Riemann integral, above equation can be rewrit- 
ten in integral form. In this way, discrete mechanics in continuous limit will give us more 
information than usual continuous mechanics [11]. For example, if we define site density at 
t — ti be p(ti) = 1/xj, above equation means the relation of site densities between time t A 
and ts is 

Here continuous canonical equations are used in integrating. It should marked that it needs 
V"p 2 + V' 2 7^ in [t A , ts], or else higher order terms in Tailer series of algorithms are needed. 
It is obvious that p(t) = constant for V(q) = aq 2 + bq, since there is V" = 0. 



V. THE ENERGY-CONSERVING MIDPOINT SCHEME AND NUMERICAL 

EXAMPLES 

In this section, we apply the energy-conserving midpoint scheme to Hamiltonian system 
that may depend on time manifestly. 

Here we consider a pendulum with periodic perturbation, whose continuous Hamilto- 
nian is H{l,6;t) = 2J ^ r2 / 2 + Mgr(l — cos#) * (1 — 0.1 * sin(0.02 * £)) where angular mo- 
mentum I and angular coordinate are the canonical variables. Take its discrete Hamil- 
tonian as ifW : H( l(k)+l 2 (k+1) , <> ik) +f +1) ■ j*±«*±i), so that energy equation (3.19) becomes 
( l (k) +l (k+i) (k) +e (k+i) _ t fc +f fc+1 N ff / ;(fc+i) + ;( fc + 2 ) g(fc+i) +g (fc+2) _ tk+tk+1 . 

n V 2 ' 2 ' 2 / V 2 ' 2 ' 2 /" 

If Mr 2 = 1 kg ■ m 2 , Mgr = 1 kg ■ m 2 / s 2 and initial values are = 0.5 kg ■ m 2 /s, 
9^ = 0.5, ti = 0s and T\ = 0.5s, then the time step-length r k := t k+ i — t k varies according 
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to the energy preserving equation (3.19). The numerical results are shown in the following 
three figures. 

Fig. 1 is a plot of time step-length r k versus time t k in < t k < 10000s. It shows there 
is a long period about ^ s that comes from perturbation —0.1 * sin(0.02 * t). 

Fig. 2 shows its detail in < t k < 350s. Besides the long period ^ s -> there is another 
short period mainly from the system without considering perturbation. 

Fig. 3 shows the detail in < t k < 16s. The upper graph shows angular coordinate 
versus time and the lower graph is versus t^. From these two graph, we find it 
gives a relative smallest time step-length when 9^> is near zero, and the relative biggest time 
step-length appears when 9^ is near an extremum. 
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FIG. 1. Time step-length versus time tk, < tk < 10000s. 
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FIG. 2. Time step-length versus time tk, < tj, < 350s. 
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FIG. 3. 8^ and r fc versus i fe , < t k < 16s. 



VI. CONCLUDING REMARKS 

In this paper, we establish an NCDC on 1- dimensional lattice with variable step- lengths. 
Based upon this NCDC, we propose the difference discrete variational principle in discrete 
mechanics for the finite time interval and apply it to both Lagrangian and Hamiltonian 
formalisms with higher order derivative with respect to discrete time. The two formalisms 
are related by discrete Legendre transformation. In order to deal the cases with higher 
derivatives, we make use of the Lagrange multiplier. Moreover we discuss some the discrete 
mechanics and symlectic algorithm in both Lagrangian and Hamiltonian formalism. We 
obtain corresponding symplectic and energy-preserving schemes and we also show that the 
necessary and sufficient conditions for the symplectic 2-form preserving in discrete mechanics 
are the corresponding Euler-Lagrange 1-form being closed [5,6]. 

When Lagrangian/Hamiltonian depend on time manifestly, we find that the discrete 
energy conservation equation in our approach is different from Lee's approach [10-12] or 
Veselov's generalized version [13]. The reason comes from whether the NCDC in base space 
as well as in fiber bundle is suitably considered. In our approach, variation problems are com- 
pletely dealt with from the viewpoint of NCDC, which is the closest analogous of continuous 
mechanics in the view of geometric-variational approach. 

Moreover, we show that the discrete mechanics in continuous limit may give some infor- 
mation of site density that cannot be obtained in usual continuous mechanics [11]. Here we 
also find that there is In p{t) ~ | \ia(y"p 2 + V' 2 ) when Hamiltonian does not depend on time 
manifestly. 

Finally, it should be mentioned that as what had been done in [5] our approach with 
the NCDC may be generalized to discrete classical field theory and the multisymplectic 
algorithm with discrete energy-momentum preserving. We will leave this subject for further 
publication. 
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